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A new computer program named CutLHCO is introduced, whose function is the implementation 
of generic data selection cuts on collider event specification files in the standardized .Ihco format. 
This software is intended to fill an open market niche for a lightweight yet flexible "consumer- 
level" alternative to the Root data analysis framework. The primary envisioned application is as 
a filter on output produced by the PGS4 and Delphes detector simulations, which are themselves 
lightweight alternatives to the Geant4 based solutions favored by the large LHC experiments. All 
process control instructions are provided via a compact and powerful card file input syntax that 
efficiently facilitates the reasonable approximation of most event selection strategies and specialized 
discovery statistics commonly employed by the CMS and ATLAS collaborations. The structure, 
function, invocation and usage of the most recent CutLHCO 2.0 program version are documented 
thoroughly, including a detailed deconstruction of several example card file specifications. The 
associated software is simultaneously being made available for free public download. 

PACS numbers: 02.70.Uu, 07.05.Kf, 29.85.Fj 



I. INTRODUCTION AND MOTIVATION 

A. Monte Carlo Event Generation 

In order to make a tangible connection between the ab- 
straction of a given theoretical construct and the actual- 
ity of detailed collider level observations, a sophisticated 
and reliable mechanism of simulation is essential. How- 
ever, the gulf separating the theoretical inception from 
the experimental inquest of a physical model can be quite 
wide. The spanning of this great divide can be logically 
partitioned into four distinct steps, each of which may be 
sequentially effected via widely established public com- 
puter code. 

The initial computational task consists of generating 
low order Feynman diagrams that may link the incom- 
ing beam to the desired range of hard scattering inter- 
mediate states. The resulting matrix elements are sub- 
sequently fed into a secondary processing phase for ap- 
propriate kinematic scaling and an implicit all-order re- 
summation of the leading order transition into batches 
of Monte Carlo simulated parton level scattering events. 
In the case that the more precise but computationally 
demanding matrix element calculation is extended to in- 
clude radiation by higher order diagrams (as is increas- 
ingly advisable outside the soft, collinear secondary emis- 
sion regime), a one-to-one matching algorithm must also 
be applied to preclude the double counting of states. Pro- 
grams popularly employed to facilitate this intricate se- 
quence of calculations include AlpGen [1| and the Mad- 
Graph/MadEvent HQ suite. 

The third tier of processing, which handles the cas- 
caded fragmentation and hadronization of the parton 



level events into final state showers of photons, lep- 
tons and mixed jets, is the domain of programs such 
as Pythia [4], Herwig [H| and Sherpa Q. Finally, 
though, a veil of obfuscation must be drawn across the 
detailed omniscience of the initial three steps, replicat- 
ing the limits on information and vulnerability to error 
of a physical detector environment. Options available 
for this fourth processing phase include the industrial 
strength Geant4 [7( simulation toolkit favored by the 
major ATLAS and CMS detector collaborations at the 
Large Hadron Collider (LHC), and also lighter- weight 
consumer-level tools such as PGS4 § and Delphes 

B. The Need for Data Selection Cuts 

The output produced by the utilities just described 
is not, however, directly suitable either for human con- 
sumption or for the efficient discernment of signal from 
background. For these purposes, a fifth processing phase 
is required - one that is sensitive to the extraordinary 
delicacy of the modern high energy experimental condi- 
tion. It is an environment in which known processes, suf- 
ficiently well understood to be relegated to subservience 
as calibration, will, by their definition as the easier tar- 
get, comprise a noise floor that tends to swamp any can- 
didate indication of new physics. The severe synchrotron 
radiation limits on light particles has moreover forced 
circular ring collider probes at the energy frontier to 
abandon the clean kinematic consumption of elemental 
electron-positron pairs for the muddled partial interac- 
tions of strongly bound quark-gluon composites. Any 
given set of final states, even assuming perfect efficiency 
in measurement, and admitting the inevitable evanes- 
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cence of the neutrino, will correspond to an innumerably 
large amalgam of unobservable internal processes. The 
stochastic variation inherent in quantum interactions will 
create false excesses and shortfalls in production which 
both mask and masquerade as the sought post Standard 
Model (SM) contributions. In this environment, discov- 
ery is only as certain as the aggregate of statistics, fastid- 
iously extracted from the tails of event distribution tails, 
accumulated with the laborious passage of time. 

The most incisive tool in the particle physicist's arse- 
nal for the clearing away of these myriad obstacles is a 
savvy application of the data selection cut. Simply put, 
it is necessary to isolate, or select, potential outcomes 
that are accessible to the desired signal, but inaccessi- 
ble to, or at least substantially unlikely for, the compet- 
ing background. It may actually be beneficial to discard 
even a large quantity of signal events from regions of 
phase space that are background dominated, in favor of 
a smaller quantity of retained events of an unusual char- 
acter that may be uniquely differentiated. Additional 
layers of selection filtering may typically be applied to 
eliminate faked signals attributable to the persistent in- 
completeness and occasional fallibility of detector mea- 
surements. Careful tailoring of the cuts to the sought 
signal may readily account for orders of magnitude of 
relative signal enhancement, and great effort is thus ex- 
pended in this pursuit, commensurate with the weight 
of potential benefit that the cuts employed may leverage 
against the great cost and effort of the project at large. 



and Delphes @, is the object of the present article. 

The CutLHCO program has evolved organically dur- 
ing extensive private use over more than a year's time for 
the rapid prototyping and application of generic selection 
cut criteria against Monte Carlo collider-detector event 
simulation data. Development under the real- world pres- 
sure exerted by a need to carefully model a substantial 
(and growing) variety of actual selection strategies from 
the ATLAS and CMS collaborations has resulted in an 
extremely flexible and powerful analysis framework that 
is capable of addressing the majority of leading contem- 
porary use cases. To cope with this extreme generality of 
function, the user interface syntax has necessarily been 
driven toward an extremely simple and intuitive form, 
with all processing instructions compactly delivered via 
a single input card file. A recently established plateau 
in stability and sophistication has presented an apt oc- 
casion for wider release of the associated package into 
the public domain under the terms of the GNU General 
Public License (l7l |. in conjunction with assignment of 
the 2.0 versioning designation. The full CutLHCO 2.0 
distribution is available for download from the author's 
personal website [HI], where it shall be updated as new 
versions become available. The main program is also 
included as an ancillary file "anc/cut_lhco.pl" with this 
document's electronic source at the arXiv.org repository. 



II. PROGRAM STRUCTURE AND FUNCTION 



C. A New Selection Cut Tool 

Selection cuts are most often implemented by the ma- 
jor detector collaborations within the Root [h]] data 
analysis framework. However, it would appear that 
there remains room within the Monte Carlo collider- 
detector simulation ecology for development of light- 
weight consumer-level event selection tools that em- 
body the analogous role with respect to Root [lOj 
that PGS4 d| and Delphes [§| play with respect to 
Geant4 0- Parallel efforts to fill this niche include 
the MadAnalysis5 (ll [package available for integration 
with the MadGraph 3fl software family, and Mathe- 
matica notebook based solutions such as the original 
Chameleon [12j package and various extensions based 
upon it. During the course of an extensive phenomeno- 
logical study with colleagues Nanopoulos, Li and Maxin 
of a particle physics model named J--SU(5) [T3 - [u^ |. it 
proved beneficial to develop and refine a proprietary soft- 
ware solution for the implementation of desired cuts, and 
the counting and compilation of the associated net statis- 
tics, in order to facilitate a clear and testable descrip- 
tion of the experimental profile that our preferred model 
might present at the LHC. A documentation of the re- 
sulting Perl program, named CutLHCO for its oper- 
ation on event files in the standardized . Ihco output for- 
mat that may be produced natively by both PGS4 [8] 



Operation of the CutLHCO program may be subdi- 
vided into three basic phases consisting of the reading of 
an event specification in the . Ihco format (|II Aj) , the fil- 
tering and reconstruction of basic physics objects ([II B|) . 
and the selection of event subsets consistent with spe- 
cific global signatures (jll Cp . An explicit deconstruction 
of the user syntax required to control these various op- 
erational phases will be provided subsequently in Sec- 
tion piip . whereas the primary purpose of the present 
discussion is a higher-level introduction to the logic, or- 
ganization and physical underpinnings of the available 
program functionality. As such, it may also serve more 
generically as a gentle contemporary review of the broad 
science and art of the collider data selection cut. 



A. The . Ihco Object Format Input 

The processing cycle of a CutLHCO program instance 
begins with the reading in of an . Ihco format event spec- 
ification file, which reasonably approximates the most 
critical aspects of object substructure that are knowable 
for a given collision within a modern high energy detector 
environment (l2| . The physical geometry of such a parti- 
cle detector generally includes an inner system of silicon 
pixel track finders with extraordinarily fine timing and 
position resolution surrounded by a network of calorime- 
ters designed to separately absorb energy through elec- 
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tromagnetic and hadronic interactions; intense magnetic 
fields curve the paths of charged particles to facilitate a 
determination of the charge-to-mass ratio, and a supple- 
mentary outer gas detector may be employed that spe- 
cializes in the identification of muons, whose calorimeter 
depositions are generally weak. An astute integration of 
the raw data data collected by the various detector sub- 
systems allows for the segregation of key particle species 
with high confidence based upon various mutual tracking 
and calorimeter signatures, and also for the efficient han- 
dling of beam pile-up conditions, multiple source interac- 
tions, intrusion from cosmic rays, and secondary vertex 
displacement (a key marker of heavy flavor physics). 

While still substantially less detailed than the corre- 
sponding Root [l(| object classification, the . Ihco format 
data subset represents a reasonable compromise be- 
tween precision and ease of use. The event constituents 
are combined into high-level physics objects consisting 
of photons, leptons (with specified flavor), and recon- 
structed hadronic jet clusters. These objects are de- 
scribed one per line in plain text, each with a variety 
of standard geometrical and kinematic parameters pro- 
vided in a reader-friendly fixed column layout. All dimen- 
sionful quantities are given in GeV-scaled natural units. 
Some degree of isolation filtering and cross-cleaning is 
typically applied prior to assembly of the high level ob- 
jects. It is likewise typical that decays in-flight of the 
tau have already been handled at a higher level of the 
analysis hierarchy, and that tau leptons reported by the 
event specification are thus only the subset deduced from 
a hadronic signature within the daughter product. While 
it is usually possible to additionally restrict the . Ihco file 
content to events passing various triggering criteria, it 
may be preferable to revert to only level-0 triggering in 
the detector simulation, so that all event selection activ- 
ity may be globally encapsulated as an exclusive domain 
of the CutLHCO program. 

In addition to the object type, each line of data in- 
cludes an additional seven pieces of information, plus 
two dummy columns reserved for future expansion of the 
format. The first four serve as a one-to-one proxy for 
the energy-momentum 4-vector of the object, translated 
into the traditional language of collider physicists. The 
pseudo-rapidity, defined as follows, is a pure function of 
the zenith angle 9, as measured from the instantaneous z 
direction of travel of the counterclockwise beam element. 

t]= -lntan(0/2) (1) 

Forward (or backward) scattering correspond to r, equals 
plus (or minus) infinity, while 77 = is an entirely trans- 
verse scattering event. The angle <f> is simply the usual 
azimuthal angle, measuring orientation around the beam 
axis. The transverse momentum P T gives the magnitude 
of the 3-vector momentum P projection that is perpen- 
dicular to the beamline. 

P T = \/p 2 + P 2 (2) 



The final parcel of kinematic data is the object invari- 
ant mass M, which is especially important for jets rep- 
resenting the composition of several lower level physical 
objects. In particular, it should be recalled that the 4- 
vector sum of individually massless objects may easily 
acquire an energy-momentum imbalance that manifests 
as a non-negligible mass-square in the invariant product. 

P^ = E 2 - P P = M 2 (3) 

The given data may readily be inverted back into a stan- 
dard 4-vector form for internal use. 

P = |p t cos0, P T sin0, P T tan ^2tan _1 e 

E= VP-P + M 2 

P^ = {E,P} (4) 

It is also useful to define a quantity AR derived from 
the pairwise kinematic descriptions of two objects that 
gives a fairly scaled measure of their relativistic "angular 
separation", in radians. 

AP = yj (At,) 2 + (Atf (5) 

The final three active data columns contain additional 
information of a non-kinematic variety, and may perform 
multiple duty by context. One column is reserved to in- 
dicate the number of identified tracks associated with the 
object. For charged leptons, the sign of the track count 
distinguishes particles from anti-particles. For hadron- 
ically reconstructed taus, the track magnitude will be 
either I or 3. For jets, it is a positive number or zero. 
The next column is used to indicate whether a given jet 
has been tagged for heavy quark flavor content (usually 
a b, but also sometimes c). For muons, this column is 
appropriated instead to provide the numeric identity of 
the jet most adjacent in AP to the muon track. It is 
often the case that this information will already have 
been employed by a muon cleaning script at a higher 
level in the processing chain to remerge poorly isolated 
muons with a parent jet cluster. If this has occurred, the 
appropriate jet kinematics will have been adjusted and 
the track count incremented for each assimilated muon; 
additionally, the integer appearing after the decimal in 
the jet track count will reflect the number of such occur- 
rences. The last data column provides the object's energy 
deposition ratio into the hadronic and electromagnetic 
calorimeters. For present purposes, it is useful to recast 
this ratio as an electromagnetic fraction (c/. Ref. [19]). 

C= (l + E^/E^y 1 (6) 

Excessively large values (£ ~ 1) may indicate that a given 
jet (especially one outside the pseudo-rapidity bounds 
|ry| > 3 where leptons and photons are typically explic- 
itly classified) is actually an electron or a photon rather 
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than a hadronic object. For muons, this column con- 
tains a decimal number, where the digits to the left of 
the decimal indicate the amount of transverse momen- 
tum carried by adjacent tracks P^ di within a AR < 0.4 
cone around the muon track, and the digits to the right 
of (and including) the decimal give a fraction represent- 
ing the transverse calorimeter energy 

£t (i.e. the en- 
ergy deposition reduced by a trigonometric factor for an- 
gular orientation) within an inclusive three-by-three cell 
array surrounding the muon track divided by the muon 
transverse momentum. Again, this information is likely 
to have already been fed into a higher level processing 
phase dedicated to muon isolation cleaning, prior to in- 
vocation of the CutLHCO script. For internal use, a 
composite transverse energy-momentum isolation ratio 
(cf. Ref. [20]) is defined as follows. 

C = (i^ + Et 1 ) /P T (7) 

Subsequent to the itemization of constituent physics 
objects, the . Ihco format specifies a final line to report 
the composite calorimeter based event missing transverse 
energy 

#t . This is, again, primarily a calorimeter based 
statistic, which is supplemented by additional kinemat- 
ics extrapolated with aid of the muon detection subsys- 
tem. In words, this quantity measures the magnitude 
of the vector sum over directed calorimeter energy de- 
posits, trigonometrically reduced to reflect only compo- 
nents transverse to the beamline. The imbalance re- 
flected by a substantial residue indicates that weakly 
interacting particles may have escaped the detector. In 
addition to the calorimeter estimate for the missing en- 
ergy vector magnitude, an angular orientation in the 
transverse plane is also provided; to be more precise, 
this orientation actually represents the angular negation 
of the described vector sum, i.e. the directionality of 
the missing energy deposition whose inclusion would re- 
close the overall sum to the expected zero magnitude. 
The topic of missing energy will be revisited in Subsec- 
tion pi CP , where the native . Ihco calorimeter based es- 
timate will be contrasted to track based estimates com- 
putable directly within CutLHCO. 

B. Object Reconstruction 

After reading the relevant content from the . Ihco input 
file, the CutLHCO program initiates an event process- 
ing phase that may be broadly labeled as "object recon- 
struction", according to controls specified by the user in 
a card file. The basic purpose of the object reconstruc- 
tion is to enforce minimum data quality characteristics 
for groups of leptons and jets, in terms of their kinemat- 
ics, geometry, proximity to other objects and multiplic- 
ity. Classified particle groups may be subsequently deref- 
erenced by a numerical identifier during the program's 
event selection phase, for testing against various discov- 
ery statistics designed to isolate new physics. Elements 



of the card file syntax that are organically encountered 
during this dialog on the object reconstruction procedure 
will be denoted by Typewriter font for future reference. 

The principal object grouping 0BJ_ALL is, as the name 
suggests, applied universally during the initial read ac- 
cess of all input physics objects. It takes only two pa- 
rameters, namely a global specification of the acceptable 
boundaries on each object's transverse momentum mag- 
nitude P T (ptm) and pseudo-rapidity magnitude \r)\ (PRM). 
Objects failing this primary filtering are permanently dis- 
carded from the event. Next, a set of instructions 0BJ_PH0 
is processed for the handling of photon type objects. In 
addition to the standard PTM and PRM filters, two addi- 
tional parameters are made available. The CUT param- 
eter specifies the acceptable range of acceptable photon 
counts after application of the previously described fil- 
ters, as a minimum and maximum; noncompliance with 
the specified limits causes the event to be rejected, or cut, 
as a whole. The JET parameter specifies how individual 
photon candidates that fail the filtering specification are 
to be handled; a value of "1" indicates that these objects 
should be reclassified as jets, whereas the default behav- 
ior "0" is simply to discard them. The initial separation 
of the charged leptons proceeds in a very similar manner, 
as controlled by individual object specifications QBJ_ELE 
(e), 0BJ_MUQ (/x), and 0BJ_TAU (r) for each generation. 
The only universal operational distinction is the addi- 
tion of an additional parameter SGN for controlling the 
lepton charge; values of and act to retain only 
the indicated sign grouping, whereas the default behavior 
"0" makes no distinction between particles and antipar- 
ticles. For muons there are two additional specifications, 
PTC and ETR, respectively describing the allowed range 
of transverse momentum P T ^ in GeV carried by tracks 
within a AR < 0.4 cone of adjacency to the muon, and 
the allowed range of dimensionless values taken by the 
transverse energy- momentum isolation ratio £ of Eq. ([7} . 

From this point forward, the outer control structures 
refer to the leptons only as a unified class. The object 
reconstruction 0BJ_LEP governs the content of the zeroth 
combined lepton classification, which constitutes also the 
default starting point of all latter groupings. Although 
the leptons are consolidated, finely grained control over 
inclusion and exclusion by flavor is retained by the ad- 
dition of a new parameter, EMT. This parameter can take 
the values "+2" or "+3" to specify the inclusion only 
of electrons, muons or taus, respectively, whereas the de- 
fault behavior "0" retains all lepton flavors. However, 
negative values "-2" or "-3" are also allowed, and 

are interpreted instead as an exclusion of the given flavor 
index; in practice, the value of "-3", which selects both 
electrons and muons for retention, is particularly useful. 
The SGN, PTM, PRM and CUT parameters carry over in the 
expected manner, but reclassification under the JET flag 
is no longer permitted. The PTC and ETR tags are also 
retained, but affect only the muonic object subset. 

The hadronic jets are treated very similarly to the uni- 
fied lepton classifications, although there are distinctions 
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in the available filtering options. The object reconstruc- 
tion 0BJ_JET governs the zeroth level jet assembly, which 
will again constitute the default starting point for all 
subsequent groupings. The PTM, PRM and CUT parame- 
ters are supplemented by four additional control specifi- 
cations unique to the jets. The parameter HFT indicates 
whether jets are required to possess a heavy flavor tag. A 
value of "1" indicates that loose tagging is sufficient, while 
a value of "2" enforces the stricter tight tagging criterion; 
the default setting "0" enforces no requirements on heavy 
flavor content. The parameter FEM specifies a range of 
acceptable values for the dimensionless electromagnetic 
fraction £ of Eq. ([6]). The tags TRK and MUO indicate the 
allowed range of counts, respectively, for individual ob- 
ject tracks associated with the jet, and poorly isolated 
muons integrated into the jet by a prior cleaning phase. 

The real power and flexibility of the CutLHCO object 
reconstruction specification comes from the ability to it- 
eratively define a practically unlimited number of addi- 
tional lepton and jet classification levels. These go by the 
control structure tags 0BJ_LEP_N and 0BJ_JET_N, where N 
is positive integer between "1" and "999". These speci- 
fications are processed in numerical sequence, and may 
refer for sourcing and comparison back to any previously 
defined specification with a lower valued index. Five new 
parameter keys are introduced, in addition to the full set 
of retained 0BJ_LEP and QBJ_JET inputs, to regulate this 
data network. The most important is the source SRC, 
which specifies the initial company of objects that are to 
be filtered. Indices provided for the lepton (jet) source 
parameter refer to previously defined lepton (jet) groups. 
Multiple index values may be listed, and the union of all 
lepton or jet objects (as the case may be) included within 
any of the referenced groups is assembled into the source 
for the current object reconstruction; negative integers 
are also allowed, and have the effect of vetoing the inclu- 
sion of all objects present within any group so referenced. 
The default value "0" adopts all objects within the zeroth 
lepton or jet classification as the current source set. A 
functionally similar parameter key CMP is available for iso- 
lating a complementary set of oppositely typed objects 
(leptons refer now to jets, and vice versa) against which a 
comparative filtering of objects may be applied. Each of 
the CMP and SRC inputs is accompanied by a corresponding 
parameter (CDR or SDR) specifying the acceptable range 
of Eq. ([5]) AR displacements between passed objects and 
either the external complementary object set or the (oth- 
erwise successful elements of the) internal source object 
set. Finally, the ANY tag allows for multiple lepton or 
jet specifications to be joined by the logical construct 
"or" rather than the default "and". The accepted in- 
put is a list of previously defined (lower valued) recon- 
struction indices of the same lepton or jet variety as the 
host, with the effect that the linked set of selections is 
passed as a whole if any single member (including the 
host) passes. This is particularly useful for constructing 
the logical negation of a compound statement, according 
to the identity not(A and B) = (not A) or (not B). In 



particular, tiered classifications may be constructed that 
retain only events failing at least one of the criteria si- 
multaneously required for admission into a prior tier. 

The final category of object reconstruction groupings 
0BJ_DIL_N facilitates the identification of dilepton object 
pairs. There is no default zeroth classification in this 
case, and the integer N again ranges between "1" and 
"999". There are four new parameter specifications par- 
ticular to the dilepton groupings, in addition to a CUT 
on the number of surviving pairs. Firstly, LEP specifies a 
single integer index corresponding to a previously defined 
source lepton grouping. Since the dileptons represent an 
independent data set, processed after all lepton classifi- 
cations are complete, the selected index does not have 
to be less than N. The zeroth lepton set "0" is an accept- 
able specification, but is not a default. Next, DLS specifies 
the target dilepton sign classification. Allowed values are 
for same-sign pairs only, "-1" for opposite-sign pairs 
only and the default value of "0" to allow all combina- 
tions. A similar specification DLF restricts dilepton flavor 
mixing, where a value of "1" allows only matching fla- 
vors to be paired, and the default value of "0" imposes 
no flavor restrictions. Lastly, the specification DMI pro- 
vides the acceptable range of the dilepton mass invariant 
M 1 ' 1 in GeV, where the invariant mass M A ' B is defined 
generically for any pair of objects A and B as follows. 

M A < B EE ^(P^ + P*) (P£+Pg) 

= \Jm 2 a + Ml + 2 (^E A E B - P A ■ P B ) 
lim => A /2|P A ||P B | (l -cosA( /3 B ' aN ) (8) 

M A =M B =0 V ' " ' \ r / 

A massless limit for the individual leptons is generally a 
fair approximation in the present case, but is not enforced 
by the program. 



C. Event Selection 

Following the reconstruction of component physics ob- 
jects into groups satisfying various data quality and 
counting requirements, the CutLHCO program pro- 
ceeds on with the second half of its analysis duties, which 
may be summarized under the heading of "event selec- 
tion". This processing phase is dedicated to isolating 
global signatures of new physics that may become appar- 
ent when the event is contemplated as a unified entity. In 
particular, various discovery statistics, either optimized 
for sensitivity to missing energy and exotic decay config- 
urations or biased against false triggers and energy mis- 
measurement, are evaluated under projection of the pre- 
viously indexed objects. Elements of the card file syntax 
that are organically encountered as a part of this dia- 
log on the event selection procedure will be denoted by 
Typewriter font for future reference. 
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Selection Tier 1 

Three elementary (but not fully independent) statistics 
are presently considered to comprise the hrst event selec- 
tion tier, each of which provides a distinct perspective on 
the basic transverse kinematics of the collision event. The 
notion of transverse variables has been broached already 
in Subsection (jll Aj) . but their ubiquity of application in 
the event selection process now comes front and center. 
The driving motivation for adoption of this perspective is 
the simple reality that the detector cannot be sensitive to 
particle flow down the beamline. Moreover, since the ini- 
tial portion of longitudinal momentum carried by the in- 
teracting "parton" constituents in a hadron collision can- 
not be ascertained, an imbalance in the detected momen- 
tum sum can only be established in the beam-transverse 
plane, where the net conserved momentum may be safely 
approximated as zero. This projection bestows an intrin- 
sic kinematic incompleteness on the transverse variables, 
by dint of which they will generically tend to indicate 
new physics via a termination point rather than a peak. 
Beyond this inherent loss of information, there are other 
obstacles that additionally complicate the direct inter- 
pretation of transverse statistics. As an important case 
in point, it is useful to consider the following definition 
of the transverse energy of a single particle. 



lim 

M=o 



M 



Pry • Prp 



E 2 - P 2 



(9) 



Firstly, note that the characteristic Lorentz transfor- 
mation response associated with energy as a time-like 
four-vector component has been corrupted through a 
quadratic reduction (while retaining the same invariant 
mass M) by the longitudinal momentum. However, this 
is not necessarily a detrimental exchange; in particular, 
the transverse energy E T has acquired by trade the in- 
teresting new property of invariance under the subset 
of Lorentz transformations corresponding to longitudi- 
nal boosts. This is not a fundamental invariance of the 
type garnered by the mass M in Eq. |3]) from contrac- 
tion under the Lorentz metric, but rather an "accidental" 
byproduct of the quiescence of P T (which lacks time- 
like and longitudinal spatial components) under orthog- 
onally directed space-time rotations; nevertheless, it is a 
very useful property that may be leveraged in interesting 
ways to link relative frames of reference along the beam- 

A B 

line. Secondly, the mutual transverse energy E T ' of two 
summed four- vector constituents A and B is not equiva- 
lent to the simple sum of individual transverse energies. 



A.B 



E' 



E L 



A.B 



P A + P L 



P, 



A,B 



P^ + P 



E 



A.B 



= t/(M 



r A.B \ 2 



(|p£' B |) 2 ?E$ + E% 



(10) 



In this sense, the transverse energy adds more like a mass 
than like an energy-momentum four-vector component. 



f T : The first, and arguably the most vital, of the di- 
mensionful transverse statistics to be modeled for event 
analysis by the CutLHCO program is the missing trans- 
verse energy $ T . The inclusive calorimeter based esti- 
mate f^f 1 of this quantity described in Subsection pi A[) 
is accessed through the selector tag EVT_CAL, and is ac- 
companied by only a single control parameter, namely an 
optional CUT in GeV. However, in addition to this intrinsi- 
cally available statistic, substantial flexibility is provided 
to the user to craft transverse momentum P T track en- 
hanced variations as follows, which may call upon 
preferred subsets of the event's physics objects. 



-1 x 



Eili ( p t cos <f>i x + Pt sin 0, y 



lim 

jy/=o 



5A,B 
T 



t,A 
T 



' jr T 



(ii) 



Keeping with the previous verbal definition of f^ 1 , ftj* 
is the magnitude of the vector negation of the sum over 
the transverse momentum P T carried by N individually 
reconstructed objects. As before, it represents the un- 
accounted contribution required in order to balance the 
open polygonal vector sum of visible transverse momenta 
to the nominally expected net value of zero. The com- 
posite missing transverse energy j£^' B adds like the com- 
posite visible transverse energy in Eq. (|10[) and is equal in 
magnitude to the mass-exclusive portion of that statistic, 
given that the missing transverse momentum vector 
is the simple inverse of the visible transverse momentum 
vector P T . To recapitulate, the merit of this statistic lies 
in the fact that key indicators of new physics, particu- 
larly the SUSY neutralino \ dark matter candidate, may 
interact only weakly and are expected to yield a missing 
energy signature upon escape from the detector. 

It is in some circles conventional to distinguish es- 
timates of this quantity that are augmented by data 
from the tracking subsystems as "missing transverse mo- 
mentum", while the term "missing transverse energy" 
is reserved for estimates based dominantly on directed 
calorimeter deposition; one may even occasionally see a 
hybrid notation that superimposes a vector symbol onto 
the slashed energy. Instead, given that there is no unam- 
biguous method for assigning a non-zero invariant mass 
lyl to any underlying physics objects that may escape the 
detector network, the presently adopted practice will be 
to designate all such vectorial space-like transverse imbal- 
ances as .pp, and all corresponding time- like three- vector 
magnitudes as $ T , resolving ambiguity as to the mea- 
surement origination by addition of the superscript "cal" 
when applicable. However, determination of a mass to be 
associated with the missing transverse energy remains a 
physically relevant issue, which is driven by the practical- 
ity that escaping particles are frequently expected to be 
produced in pairs that may be entirely kinematically de- 
coupled from each other by intervening steps in the decay 
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cascade, implying that a substantial portion of the un- 
seen tracks associated with the true missing momentum 
vector may internally cancel in the observable sum. In 
fact, this dilemma forms the motivation for certain of the 
specialized discovery statistics to be introduced shortly 
as the third event selection tier. 

The zeroth incarnation of the missing transverse en- 
ergy f T to be computed internally is labeled for input 
purposes as EVT_MET, and includes the tracks of all con- 
stituents of the event that survive the primary 0BJ_ALL 
classification. It is made available for reference to latter 
selectors by default, with no need for specific invocation 
by the user. Again, an optional CUT may be specified 
in GeV. Additional indexed EVT_MET_N selectors, with N 
an integer in the range "1" to "999", are provided for 
evaluation of the missing transverse energy carried by 
various object subsets. Two additional parameter spec- 
ifications, LEP and JET, allow sourcing from previously 
defined 0BJ_LEP_N and 0BJ_JET_N object reconstructions. 
At least one of these inputs must be set equal to a single 
integer back-reference N; a value of "0" is allowed, but is 
not a default. The missing energy event statistics may 
themselves be made available for back-reference by their 
integer index to the MET parameter of subsequent selector 
definitions, with EVT_MET taking the implicit valuation of 
N = "0". As a special case, the calorimeter based estimate 
EVT_CAL is also automatically made available, under the 
index assignment N = The directional characteris- 

tics of each computed missing transverse energy statistic 
are retained internally for reference by future selectors. 

H T : A second event selector of great importance is 
the scalar sum on transverse energy H T , which is defined 
as follows and indicated programmatically as EVT_MHT. 

lim => Ei-x Pi (12) 

A massless limit is imposed by default, but may be over- 
ridden by setting the MAS parameter to "1". The key dis- 
tinction between H T and the missing transverse energy 
(in the massless limit) is that the absolute value is ap- 
plied prior to the summation rather than after, so that 
vector cancellation is avoided. As such, this statistic is 
instead representative of the net energy carried by the 
visible event constituents. In fact, it represents precisely 
the massive transverse longitudinal scalar possessing the 
property of linear composition that eludes the transverse 
energy as defined by Eq. (fTUI) . 

Hrr*'^ = + Hr£ (13) 

In this sense, the massive limit might be considered some- 
what more deeply consistent, insomuch as the desired 
associative property may be disrupted by compound al- 
ternation of mass-zeroing with four-vector summation. 
As before, H T is automatically computed for all primary 
event objects in the zeroth incarnation, and an optional 
CUT is allowed on the result in GeV. Likewise, indexed 



EVT_MHT_N H T selectors are available as before for treat- 
ing LEP and JET object subsets. 

Mt : The missing and visible energy magnitudes are 
combined in a third selector called the transverse effective 
mass Af*, which is defined as the following simple sum, 
and indicated programatically as EVT_MEF. 

Mf = f T + H T (14) 

Once again, is automatically computed for all pri- 

mary event objects, with a single optional parameter CUT. 
Additional indexed EVT_MEF_N M£ selectors are likewise 
provided as before, but the supplemental parameter spec- 
ifications are different in this case. Since M^ S will be con- 
sidered as a dependent construction based upon f] T and 
H T , the available inputs are instead MET and MHT, each of 
which must specify an integer back-reference N to a pre- 
viously defined incarnation of EVT_MET_N or EVT_MHT_N; a 
value of "0" is allowed in either case, but is not a default. 

Selection Tier 2 

A second tier of derivative event selectors is consti- 
tuted from various comparative combinations of the three 
previously defined statistics. No default analysis is per- 
formed in these cases, and an index N in the usual "1" 
to "999" range is always explicitly specified. The first 
such statistic, EVT_RET_N, consists of the dimensionless 
comparison 

(#t/#t) of two missing energy magnitudes. 
This quantity might be used, for example, to protect 
against the danger that reinclusion of discarded soft or 
high pseudo-rapidity jets could substantively alter the 
observed momentum imbalance. Back-references to a 
pair of previously defined missing energy indices must be 
provided to the inputs NUM and DEN, for use in the numera- 
tor and denominator, respectively. The CUT specification 
accepts pure numerical values for evaluation against the 
computed ratio. A similar pair of selectors EVT_RHT_N 
and EVT_REF_N are provided for analysis of the ratios 

{$t/V[Ht]) an d (^TA^T ffD )> which instead reference 
the scalar momentum sum and effective mass in the de- 
nominator, respectively. The square root employed in the 
EVT_RHT_N example evinces the expected Gauss/Poisson 
fluctuation size, so that the statistic provides a metric for 
the relative significance of the missing energy magnitude. 
The former NUM, DEN and CUT inputs extend in a straight- 
forward manner to each of these latter two cases, al- 

l li 

though one should take note of the curious GeV ' units 
for the EVT_RHT_N selector output. The final statistic to be 
defined, EVT_DET_N, represents the absolute vector differ- 
ence — frl between two distinct population studies 
(A,B) of the missing transverse momentum. Its function, 
once again, is to protect against false indications of es- 
caped particles. In addition to a CUT specification in GeV, 
two parameters, ONE and TWO, are provided for specifying 
(in either order) a pair of back-references to previously 
defined missing energy indices in the usual manner. 



s 



Selection Tier 3 

The third tier of provided event selectors represents 
a sampling of somewhat specialized statistics that are 
in contemporary use by the ATLAS and CMS collabo- 
rations. As with the second tier, no default computa- 
tion is performed in these cases, and an index N in the 
usual "1" to "999" range is always to be explicitly speci- 
fied. A robust variety of analyses are provided with the 
presently documented CutLHCO 2.0 release, although 
the creation of additional subroutines to supplement the 
generality of this coverage represents a potentially active 
source of ongoing program development. The author wel- 
comes useful suggestions for expansion in this regard, and 
encourages users who write their own program extensions 
to resubmit them for inclusion into the main program 
trunk. 

M%- , : The first specialized statistic to be discussed is 
the transverse mass computed for a leading lepton 

object reconstruction and the missing transverse momen- 
tum, which shall be denoted for programmatic purposes 
as EVT_LTM_N. This statistic has been employed frequently 
by the ATLAS collaboration (2l| - [23j . althou gh e ssentially 
similar channels are also of interest to CMS [24j . For any 
two objects A and B, the transverse mass statistic is de- 
fined as follows 1 25]. 
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(15) 



Note that the referenced transverse energy E T of each 
constituent object is defined in a recursively consistent 
manner according to the prescription by Eq. ©. Also, 
the simple sum of individual transverse energies that is 
featured in the Eq. (|15[) definition mimics the £f A ' B pre- 
scription from Eq. (fT5|) rather than the behavior of E^' B 
from Eq. (jTDJ) . In fact, by inspection, the mismatch be- 
tween these quantities is precisely also the mismatch be- 
tween the pairwise invariant mass M A ' B of Eq. §5§ and 
its transverse variant M, 



A.B 



The Af T transverse mass distribution will have an 
cndpoint corresponding to the mass of the presumed 
common parent of the included lepton and missing 
energy candidate, as may be verified by noting that 



A,B 



A B 

T < M ' . This potentially facilitates use as a tool 



Mr 

for the detailed mass measurement of SUSY particles, al- 
though the key present utilization is simply as a discov- 
ery statistic, where a lower mass bound in the vicinity of 

A B 

100 GeV is typically imposed. The transverse mass 
is distinguished as a Lorentz invariant under the subset of 
transformations consisting of longitudinal boosts for the 
same trivial reasons that the definition of E T from Eq. ^ 
was itself so recognized. An individually massless limit 



for both the lepton and missing transverse momentum is 
enforced in the present program implementation. 

The available control parameters include the usual CUT 
selector in GeV units, and index back-references for the 
included LEP object set (from which the leading lepton is 
extracted) and the counterpart MET vector. 

M^: A generalization of the previously described 
transverse mass statistic referred to as Mp| or the jet "s- 
transverse mass", and which will be addressed program- 
matically as EVT_JSM_N, has been adopted by the CMS 
collaboration [26| in order to better treat the SUSY sce- 
nario where pairwise sparticle production leads to dual 
decay chains that each manifest a missing momentum 
signature. First introduced in Ref. [27|, the s-transverse 
mass construction is based upon the observation that 
each such visible shower may be convolved with its af- 
filiated contribution to the missing momentum to yield a 
transverse mass, as defined in Eq. (fTS"|) . that is bounded 
above by the mass of the system's parent particle. How- 
ever, since only the unified missing transverse momen- 
tum vector fj, is experimentally available, further nu- 
merical analysis requires the specialization to a suitably 
optimized location in the distribution of this sum be- 
tween the pair of missing momentum candidates. The 
M T2 prescription is to assume an equivalent structure 
for each decay chain, so that the larger of the two trans- 
verse masses may be used imply a lower bound on a 
common parent species, subsequent to minimization of 
that maximum value with respect to all consistent 
partitions. If event pollution by upstream decays of the 
initial state prior to the targeted pair production can be 
neglected, the minimization may actually be performed 
analytically [28[. If the mass of the escaping particle 
species x IS additionally set equal to zero, the following 
simple expression is realized. 



A'p = E^ E^ 



pA p-B 



M, 



A,B _ 
T2 — 



lim 

M A =M R =0 



2|PtI 



P B \ 



1 + cos A(j> 



B,A 



(16) 



Relative to Eq. ([15]) , one immediately notices the curious 
emergence of a Euclidean signature for the inner prod- 
uct. This difference may be traced to the fact that M T2 
employs two visible systems rather than a single visible 
system in conjunction with associated missing energy; es- 
caping missing energy is expected to be antiparallel to its 
conjugate visible shower, inducing a phase shift by it ra- 
dians that effectively reverses the sign of the trigonomet- 
ric factor. In order to apply the formula from Eq. (|16p . 
the reconstructed objects taken to constitute the dual 
visible systems must be specified; following the CMS 
lead [26j, CutLHCO computes the s-transverse mass 
M^'l of a pair of optimally separated pseudo-jets in the 
massless limit. The "event hemisphere reconstruction" 
procedure (c/. Ref. [2{J, Section 13.4) begins by establish- 
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ing a pair of pseudo-jet axis "seeds" consisting of the two 
individual massless jets that possess the largest relative 
4-dimensional invariant mass. These axes Pf!~ are itera- 
tively reconstructed as the 4-vector sum (with rezeroed 
mass) of all massless jets that are assigned to them in 
a given round by merit of the smaller Lund-inspired [4| 
dimensionful distance measure (£aM/vb)/(-^a + E B ) 2 . 
The cycle terminates upon the retracing of a previously 
encountered hemisphere partitioning, which may (but 
does not necessarily) indicate convergence to a stable 
fixed point assignment of the pseudo-jet axes. Although 
originally introduced as tool for deducing the mass scale 
of decaying SUSY particles, the s-transverse mass (like 
the transverse mass itself) is currently enjoying service as 
a discovery statistic, with variously tuned cutoff thresh- 
olds ranging from around 100 GeV up to about half a 
TeV [26j. The available control parameters include only 
a CUT selector in GeV units and an index back-reference 
to the included JET object set. 

AE^ U : The next specialized discovery statistic to be 
made available is the CMS jet and dilepton-Z transverse 
energy balance AE^ U [3(3|, which will be referenced by 
the tag EVT_JZB_N. Searches implementing this statistic 
are designed to target events featuring a leptonically de- 
caying Z boson, which may be inferred by the signature 
of residual opposite sign electron or muon pairs. The pro- 
duction of Z bosons, in conjunction with jets and missing 
energy, is a potentially regular byproduct of SUSY tran- 
sitions from a heavy to light neutralino mass state. More- 
over, this selection naturally suppresses large portions of 
the SM event background, including the difficult QCD 
multijet component. Troublesome opposite-sign dilepton 
background competition persists from the dual leptonic 
decay of pair produced top quarks tt, which may be ob- 
viated by leveraging the absence of flavor correlation in 
this process. The AE^ statistic constitutes an active 
discriminant against interference from SM Z bosons plus 
initial state jet radiation by highlighting the kinematic 
correlation that may persist between the jet track mo- 
mentum imbalance and the observed dilepton pair, if they 
indeed were borne of a common antecedent. It is formally 
defined as follows, employing a massless limit for each of 
the composite iV-jet and dilepton objects. 

ae% u ^ lEjLi^l-IEL^I 

« IfT + ^Tl-I^Tl (17) 

The second line follows as a precise equality if the full 
missing momentum signal is attributable only to the 
dilepton and included jets. The unwanted background 
events are expected to be symmetrically balanced about 
the zero line of this statistic, with random fluctuations 
driven in either direction by faked missing energy. For the 
targeted SUSY events, directional collinearity between 
the Z boson and the escaping light neutralino will sys- 
tematically bias AE^ U toward a positive value. The 
active signal region is typically restricted to that sub- 
set possessing a dilepton invariant mass M 1 ' 1 (cf. Eq. [5]) 



within about 20 GeV of M z . Additionally, one generally 
requires the presence of at least three jets, to ensure de- 
coupling between the orientation of potential momentum 
mismeasurements and the Z boson track. Along with a 
standard GeV CUT selector, the user must specify back- 
reference indices for a pair of previously defined JET and 
DIL object sets. 

& a": A trio of SUSY search statistics cur- 
rently employed by the CMS collaboration [jy], H2] are 
the razor variables Af TR and M R , and their dimension- 
less squared ratio, denoted here as a R . Based upon a 
suggestion originally proffered in Ref. j33[, the pair of 
massive razor quantities represent independent estima- 
tors, the former transverse and the latter longitudinally 
inclusive, of the scale (m~ — m|)/m^ underlying com- 
pound SUSY event chains, where is the common mass 
of an initial hard squark or gluino pair production event 
and mj represents a pair of weakly interacting particles 
from the dual decay cascade that escape the detector un- 
seen. In practice, the jet razor mass variable M R , whose 
modeling is invoked by the tag EVT_JRM_N, may be availed 
to isolate the scale of prospective new physics in conjunc- 
tion with application of the a R ratio, which is provided 
the independent identification tag EVT_ALR_N, as a "ra- 
zor" discriminant against SM background competition, 
primarily from QCD di-jets. The M TR variable, which is 
dependent in this point of view, is not given its own ac- 
cessor tag, although the corresponding numerical value 
may be easily reconstructed as needed from the other 
two. The razor statistics confront the same innate diffi- 
culty as the s-transverse mass analysis, namely that 
the history of an event with a dual x detection failure is 
essentially under-constrained by the residual visible par- 
ticle content, generating fundamental ambiguity within 
any reconstruction effort. The strategy for surmount- 
ing this obstacle that is embodied by the razor analysis 
rest upon an approximation of the opaque center of mass 
and pairwise q rest frames (which will be identical in the 
limit of pair-production at the kinematic threshold) via 
the expression of a transparent i-axis boost into the "ra- 
zor frame", wherein the visible center of mass energy may 
be estimated in terms of the explicit longitudinal invari- 
ant M R composed from observations projected solely out 
of the laboratory frame; the garnered longitudinal invari- 
ance protects the integrity of the M R distribution peak 
even if the boost into the true center of mass frame is non- 
trivial. Construction of the razor frame is initiated by 
partitioning all relevant physics objects into two massless 
pseudo-jets A and B representing the observable remnant 
of the putative SUSY pair production event. The histori- 
cally favored methodology is replicated by election of the 
object partition that minimizes the sum of pseudo-jet in- 
variant mass-squares. The jets are taken to be individu- 
ally massless and, following their reassembly, the masses 
of each composite pseudo-jet are likewise zeroed by an 
appropriate compensation of the energy. The specific re- 
quired longitudinal boost is then guaranteed to exist and 
is expressible as (3 = (P^ + Pf)/{E A + E B ), assuming 
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vanishing of the transverse initial state radiation com- 
ponent. The razor variables themselves are defined as 
follows. 



M, 



A.B 
TR 



{if tI (I-PtI + I-PtI) - f t • (Pt + Pt)}/2 
W(< A ) 2 + (M T ^ B ) 2 



lim 

jy*=M A =M B =0 



A.B 



(E A + E B f - (P A + Pff 



A,B \2 



"R 



(M 



(M TR /M R )^ 



E! 



(18) 



In the transverse razor mass, it should be noted that 
a legitimate missing momentum, even if constructed in- 
clusively from the event at large, should be very nearly 
equal and opposite to the pseudo-jet vector sum, making 
the second term under the radical positive, and equal in 
magnitude to the missing energy squared. If the pseudo 
jets are close to parallel, then each term in the sum will 
moreover be equal, and will take the scale of the 

missing energy as a whole. Conversely, if the pseudo-jets 
are back-to-back, as is increasingly likely for QCD di- 
jets, a partial cancellation will be exhibited, even in the 
presence of faked missing energy. If there is no missing 
energy signal at all, then M TR will be identically zero. 
The first term under the radical may be recognized 
as the Eq. ([5]) invariant mass-square of the merged visi- 
ble pseudo-jet system, and the expression as a whole is 
identical to the definition of the transverse energy, as in 
Eqs. (|9ll0p . for this same system. Since both mass es- 
timators are designed to highlight the same underlying 
physical scale, the ratio should be expected to peak 
at order one for signal events; specifically accounting for 
"geometric" factors in the respective Eq. (fT5)) definitions, 
the fact that M TR (like all transverse statistics) identifies 
a kinematic endpoint, and a small relative Lorentz trans- 
formation factor 7 > 1, the signal peak will tend to occur 
closer to a value of 1/4. The SM QCD distribution will 
indeed peak at zero, allowing for a clear-cut excision of 
the encroaching background competition. The available 
control parameters for each of the individually address- 
able statistics includes a CUT selector in GeV units and an 
index back-reference to the included LEP and JET object 
set; for EVT_ALR_N the user may additionally include a 
reference MET to the preferred missing transverse momen- 
tum, which will instead be rebuilt from the local lepton 
and jet content if an index is not provided. 

c4 J : A replication of the CMS a T ratio QjJlilll 
is provided for use under the identification handle 
EVT_ALT_N. This quantity was devised [HI to help dis- 
tinguish actual missing transverse energy from detector 
mismeasurement. It is formally defined as follows, where 

B A 

AH T ' is the minimal positive difference that may be 
realized between the pair of _ff T sums corresponding to a 
grouping of all included objects into two effective pseudo- 
jets A and B, when considering all such possible parti- 
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If there is no mismeasurement or true missing energy, 
the value of a^ J will just be 1/2. For energy mismeasure- 
ments of otherwise anti-parallel pseudo-jet pairs, subtrac- 
tion of the nonvanishing scalar difference AH^ J will tend 
to drive ai^ 3 below the midline. In this case AH T = f! T , 
but the squaring of the small (typically <S 1) factor in 
the denominator renders it less significant. Conversely, 
genuine missing energy, as manifest in the departure 
from pseudo-jet anti-parallelism, will imbalance the vec- 
tor sum within f T more so than the simple magnitude 
difference AH T , tending to create a contrasting elevation 
in cxij? above one-half. In addition to a dimensionless CUT 
selector, the user may specify index back-references for 
the included JET object set, and the referenced MET and 
MHT values. The latter two quantities will be recomputed 
from the jet content itself if omitted. Additionally, the 
same MAS tag available to MHT itself is also available here 
to indicate a preference for the mass inclusive determina- 
tion of H T . As an editorial aside, it has been observed by 
the author's research collaboration [T^ . [lij that the a T 
statistic may become unsuitable in the soft, high mul- 
tiplicity jet regime, as the large array of finely grained 
pseudo-jet permutations makes quite likely that a bal- 
anced scalar sum AH T ~ might be achieved, creating 
an overall suppression of the final computed value. 

A</> J; ^: The next event selector to be discussed is a 
relatively simple angular difference statistic A(f>j ^ em- 
ployed by the ATLAS H3 collaboration. It will be la- 
beled EVT_MDP_N, and is defined as follows. 
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To summarize, A<j) 3 ■ ^ simply isolates the minimal angu- 
lar separation in the transverse plane between a set of 
N reconstructed objects and a fixed missing energy vec- 
tor. It is designed to protect against spurious indications 
of missing energy associated with detector mismeasure- 
ment. The user may specify a CUT on the resulting to 
7r radian angle, as well as index back-references for the 
active JET object set and the static MET vector, the latter 
of which will be reconstructed from the included jets if 
omitted. To reduce the likelihood of randomly generating 
a dispositively small value for Acftj* concurrently with a 
legitimate missing energy signal, this statistic is typically 
computed for no more than the four hardest jets. 

A</>*£: The final event selector to be modeled is a 
CMS statistic [l9[ that is quite similar to the just dis- 
cussed angular difference A^- *, but with a somewhat 
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more complicated implementation. This "biased" az- 
imuthal difference A<^>,- is formally defined as follows, 
and will be associated with the program invocation tag 
EVT_BDP_N. 
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In words, A0j # is the minimal value of the set of N an- 
gular separations in the transverse plane between each 
reconstructed object's own momentum vector and the 
missing transverse momentum composed from all other 
objects. It is expressed as an angle in radians, in the 
range to tt. Intuitively, it registers whether the elonga- 
tion of a single momentum vector is capable of substan- 
tially repairing an imbalance in the missing transverse 
momentum. If a single jet mismeasurement is indeed 
dominantly responsible for a false missing energy signal, 
then A(f>j * should register close to zero. In addition to 
an angular CUT selector in radians, the user may specify 
index back-references for the included JET object set, and 
the baseline MET vector. If the latter quantity is omitted, 
it will be recomputed from the jet content itself. This 
statistic may again suffer in the case of very high jet 
multiplicities; given a wide selection of randomly oriented 
jets, it becomes quite likely that the angular orientation 
of at least one such member might be sufficiently well 
azimuthally aligned with the true missing energy track 
that its rescaling could apparently rebalance the event. 



III. PROGRAM INVOCATION AND USAGE 
A. Card File Command Syntax 

All object reconstruction and event selection opera- 
tions available to the CutLHCO program are controlled 
by command lines specified in an external card file. An 
overview of the logic and scope of the available function- 
ality was provided previously, in Section ([11]) . Five exem- 
plar cards from the dozens of working selection strategies 
accompanying the current program distribution [18| are 
documented in the forthcoming Section pVI) , with a line- 
by-line textual deconstruction of the specific command 
formatting and intent provided in each case. A compact 
yet comprehensive synopsis of the command lexicon is 
provided in the Appendix. The purpose of the present 
section is a description of the uniform aspects of the card 
file command syntax and grammar. 

Each processing instruction is placed on its own line in 
the card file and begins by invoking the name of an object 
reconstruction (0BJ_ABC_N) or event selection (EVT_ABC_N) 
identifier, where ABC represents a three character alpha- 
betic code unique to each command family, and N is an 
integer within the range "0" to "999". Omission of the 
integer N and the associated trailing underscore is func- 
tionally equivalent to a selection of 0. This zeroth integer 
specification is generally reserved for quantities that are 
automatically evaluated by the program, irrespective of 



user input, although it may still be possible to modify the 
default parameter configuration in these cases through 
an explicit reference. However, not all input tags pro- 
vide a meaningful zeroth configuration, and, conversely, 
not all of those tag specifications that do automatically 
provide for a zeroth configuration additionally support 
multiple numerical versioning. Any lines of input that 
do not begin in the manner described are disregarded as 
comments. Likewise, the inclusion of a "#" hash symbol 
at the end of a valid input line is interpreted to initiate 
a comment enveloping all trailing content. The order in 
which lines are specified in the card file is not material to 
the sequence in which they are ultimately evaluated; the 
order of program execution is instead established by the 
fixed operations sequence described in Section ([II)). and 
within that framework, whenever there might otherwise 
be ambiguities, by an ascending numerical sort on the 
associated numerical tags. 

The parameter specification input for each command 
line is separated from the leading identifier name by an 
equals sign "=". This parameter specification consists of 
a list of key and value pairs, individually joined by a 
colon " :", and separated from adjacent pairs by a comma 
Each parameter key consists of a three character 
alphabetic string that uniquely specifies the role of the 
subsequently input value. The valid parameter specifi- 
cations vary according to the specific object reconstruc- 
tion or event selector context of each command identi- 
fier, as detailed in the document Appendix. The order in 
which a given command line's parameter specifications 
are provided is a matter of the user's discretion and has 
no impact on the logical flow of the selection algorithms. 
The value assigned to each key may consist of an integer 
or floating point number (including signed values), the 
character string "UNDEF" (to specify an undefined or null 
input), or a list consisting of elements drawn from the 
prior classes, individually separated by commas and 
enclosed in square brackets "[ ... ]". The set of value 
inputs that may be considered valid or useful is, again, a 
function of the context provided by the active parameter 
key. Unspecified values default to "UNDEF", or if the ex- 
pected input is an integer to be chosen from a predefined 
list of options, the value "0". Parameters designed to ac- 
cept multiple inputs will treat an isolated value as the 
leading list element, with all other terms undefined; like- 
wise, there is no penalty for enclosing a single expected 
input in square brackets, creating a list composed of a 
single entry. 

A frequently employed value formatting idiom is the 
"[Min, Max]" pair, used to indicate a numeric range. So 
long as Min < Max, the set of matched numbers sim- 
ply consists of those values inclusively bounded by the 
specified range. An undefined value for Min is treated 
as indefinitely small, and an undefined value for Max is 
treated as indefinitely large. If no numerically valid lim- 
its are provided, then all values match. If the upper 
and lower boundaries are numerically equivalent, then 
only that single value matches. However, a subtle ad- 
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ditional functionality is accessed under the circumstance 
that Min > Max. In this numerical match is achieved 

if the comparison value is either at least as large as Min or 
at least as small as Max; in other words, the interval from 
Max up to Min is rejected, exclusive of the boundaries. Of- 
ten, an additional instruction flag is accepted from the 
third position of the described value list, which may be 
used to indicate various processing preferences, accord- 
ing to context. The availability and function of this third 
input, whose uses include sorting objects that match 
the specified range criteria, terminating object match- 
ing upon an initial failure, and clipping the count of ob- 
jects satisfying some independent selection, is efficiently 
summarized as part of the encompassing command spec- 
ification reference provided as a document Appendix. 



B. Running the Program 

The CutLHCO "executable" is delivered as a single 
Perl script named cut_lhco.pl; to be more precise, since 
Perl is an "interpreted" language, this file serves both 
as the program source document and runtime portal. A 
benefit of this paradigm is inherent platform indepen- 
dence, making the CutLHCO program ready for im- 
mediate use, without the need for installation or compi- 
lation, on any computer with a reasonably modern perl 
environment; this should automatically include almost 
all systems in current service that run a Unix/Linux fla- 
vor variant, including Mac OS X, and free installations 
for Windows are also readily available on the web. The 
detailed instructions to follow will assume a Unix styled 
directory structure and command line shell. 

The CutLHCO distribution package is available for 
download from the web [lH, and is delivered as a tar- 
gzipped file bundle. To unpack the distribution, the user 
should run the commands "gunzip cut_lhco.tar . gz" and 
"tar -xf cut_lhco.tar" in sequence from the command 
line. The directory from which CutLHCO is called must 
possess two specific subdirectories named Events/ and 
Cards/ ', which will house, respectively, the . Ihco format 
input files and the card files controlling the application 
of selection cuts. These directories are created automat- 
ically upon unpacking the program distribution as de- 
scribed. Users of the MadGraph |3| simulation envi- 
ronment will recognize this directory structure, which is 
in fact appropriated to streamline integration with that 
framework. At the user's discretion, the main program 
file may likewise be grouped with any existing scripts 
and executables in a common bin/ subdirectory. Exe- 
cution of the script can then be as simple as entering 
" ./cut_lhco.pl", or "./bin/cut_lhco.pl", at the termi- 
nal window, as the case may be; however, there are 
several optional command-line parameters that the user 
may wish to subsequently specify, which shall now be de- 
scribed in the sequence of their usage. Any such instruc- 
tions are typed directly after invocation of the program 
name, separated from it and from each other by an empty 



space; any parameter position may be effectively skipped 
by providing the value "UNDEF". 

The initial parameter input is associated with the file- 
name that one wishes to process, and it takes a de- 
fault value of 11 data" . Any processed output will ul- 
timately be placed into a corresponding file named 
11 Events /filename, cut". The referenced input file must 
physically exist within the Events/ directory, either as 
" filename. Ihco" or as" filename _pgs_events. Ihco" . However, 
the latter scenario actually allows for two possible points 
of flexibility. A group of files within Events/ that fol- 
lows the specified root name with an extra underscore 
and a positive integer, as may be produced by the Mad- 
Graph [3j helper script "multi_run" , will be treated as 
a family of batch jobs to be merged for analysis. Also, 
any matching files with a trailing extension of .gz will 
be automatically unpacked with the gunzip utility prior 
to subsequent processing. Additionally, the more simply 
named "filename. Ihco" file will be created, and filled with 
the appropriate merged content. The next parameter 
input position is associated with the cardfile governing 
event selection, and has a default value of "cut_card". 
The referenced file must likewise exist within the Cards/ 
directory, as "cardfile.dat". 

The third parameter position allows the user to spec- 
ify a numerical cross section in pb for the physical pro- 
cesses under consideration. Specifying the negative cross- 
section value causes the program to search for a file 
named "filename pythia.log" (or the set of integrally la- 
beled file variants) in the Events/ directory, from which 
the cross section will be dynamically extracted, along 
with a statistic for the fragmentation cut percentage. An 
unrecoverable error signal is issued if the expected file 
(or files) are not found. If multiple records are queried, 
the cross section and fragmentation cut values are av- 
eraged, with equal weight assigned to each individual 
file. The fourth position designates a target luminos- 
ity in pb" 1 to which the surviving event count should be 
scaled. Since this calculation is dependent upon a known 
cross section, the reading of this input is skipped unless 
the third parameter was defined and non-zero; if required 
but not defined, the luminosity target is set to a default of 
1,000 pb" 1 ee 1.0 fb _1 . The fifth and final input, which 
may actually shift into the fourth position as described, 
represents an integer variable for controlling the output 
of detailed event-by-event statistics. The default value 
"0" suppresses this output, while "1" enables it, and "2" 
additionally implements a resorting of the original event 
sequence by the numerical value of each data column, 
from left to right; this sorting process may take an ap- 
preciable amount of time for large quantities of surviving 
events with large multiplicities of statistics. 

To provide a concrete initiation example, the com- 
mand "./cut _lhco.pl ttbar cut_card_LHC 79.8 4500 1" 
will, given placement of the main executable within the 
working directory, cause a file named "ttbar. Ihco" (or its 
permissible variants) within the Events/ directory to be 
processed according to the selection cut strategy de- 
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scribed in the file " Cards /cut card _LHC.dat" , assuming a 
process cross-section of 79.8 pb, with results scaled up to 
a sample size of 4.5 fb -1 , and all event statistics reported 
but not sorted. 



C. Program Output and Interpretation 

The output . cut file described in the prior section pro- 
vides a summary report of the CutLHCO runtime en- 
vironment and any actionable statistics generated from 
the surviving events. It opens by stating the number of 
. Ihco files merged during the process initiation (if there 
were multiple inputs) , the total number ./V of events pro- 
cessed, the percentage of events failing fragmentation 
cuts (if harvested from the Pythia [4| log), the source 
event cross section a in pb and associated simulation lu- 
minosity C = N/cr in pb" 1 (if available), and the target 
luminosity and associated event rescaling (if non-trivial) . 
These particulars of the operational context are followed 
by a presentation of key results and analysis, beginning 
with the all-important scaled count of events surviving 
the selection cut regimen. Next, the percentage of events 
cut during the selection is tabulated, as a unified total 
and also as a statement of activity for each individually 
defined cut; the latter set of statistics is supplemented 
by the percentage tally of circumstances in which each 
specified cut constitutes a uniquely enforced rationale for 
event exclusion. Each object or event command provid- 
ing a value for the "CUT" selector that is not function- 
ally equivalent to "UNDEF" will be included in this report, 
referenced by its unique alphabetic identification code 
(omitting the "DBJ_" or "EVT_" prefix) and its numerical 
sequencing tag. 

If a detailed itemization of results was specified during 
the program initiation, then each surviving event will be 
given a line in the output that inventories its response to 
each defined selection cut query. This is particularly use- 
ful for tabulating histogram reports in various discovery 
statistics, or for partitioning the counts of events satisfy- 
ing subordinate selection conditions. It may also be use- 
fully applied to optimization of the underlying selection 
strategy itself, facilitating the identification of inflections 
or knees in the event distribution, potential advantages 
over (or novel tools for the suppression of) background, 
and selection redundancies or inefficiencies. To include 
a given output statistic in this report without actually 
imposing a cut, a trivial value should be provided for the 
CUT selector that is still technically defined, e.g. an inclu- 
sive lower bound of "0" for an inherently positive quantity 
such as P T . A helper script histogram.pl is provided in the 
main directory of the CutLHCO distribution [l8| that 
is well suited for this analysis task. This script concur- 
rently processes all . cut format files within the Events/ 
subdirectory, and place the results into a file called " his- 
togram. txt" in the working directory. 

The histogramming script accepts up to six optional 
space-separated parameter inputs after the primary in- 



vocation ". /hi st ogram.pl", and again permits the usage 
of "UNDEF" for reversion to system defaults. The first input 
should be the name of the event detail column heading 
from the relevant . cut files that the histogram should 
be based upon. The default behavior is a simple single- 
column histogram on the net event count; although this 
repeats raw information already provided in the under- 
lying report, it may still be useful for rapidly assembling 
the counts of multiple event samples into a unified pre- 
sentation. The second input is a numerical specification 
for the binning width; the default value of "0" indicates 
a single, inclusive bin. The next two inputs numerically 
specify the minimum and maximum cutoff values for in- 
clusion in the histogram; the lower boundary defaults to 
"0", and the upper boundary defaults to "UNDEF", which 
is interpreted to be indefinitely large. The fifth input is 
a integer list selector for the numerical label that should 
be assigned to each bin; "0" (which is the default) la- 
bels each bin by its left-most (minimum, inclusive) value, 
while "1" substitutes reporting of the central bin value, 
and "2" instead calls for the right- most (maximum, ex- 
cluded) boundary value. The last parameter input is a 
binary flag for governing the histogram output format; 
the default value of "0" produces standard tab-delimited 
text that is suitable for passing on to spreadsheet based 
statistics packages, while the alternate specification of "1" 
will yield a file employing the array input syntax favored 
by Mathematica. 

Additional helper scripts for interpreting the statistics 
of surviving groups of events and for automating the gen- 
eration of supervisory shell scripts to shepherd the pro- 
cessing of large batches of parallel analyses are available 
to interested parties upon personal request. No public 
documentation currently exists for these tools, although 
the curious are welcome to explore and modify any pro- 
vided code for adaptation to their own purposes, under 
the same terms [ljj as the rest of the package. In general, 
the author is quite open to private inquiries, suggestions, 
and requests for support by potential CutLHCO users 
on any topic regarding application of this software to the 
physics analysis of collider data selection cuts. 



IV. CARD FILE CASE STUDIES 

The most efficient way to communicate a working vo- 
cabulary of CutLHCO commands to the potential user 
is probably by direct example. The general syntax tu- 
torial from Section (|III[) . along with the summary of 
available parameters in the document Appendix, are thus 
supplemented in the present section by five specific case 
studies of card files designed for the replication of actual 
LHC SUSY searches, three from the CMS collaboration, 
and two from ATLAS. The coded instruction sets are are 
transcribed in Typewriter font within individually boxed 
panels, and a line-by-line deconstruction of the event se- 
lection intent accompanies each such inset. A growing 
library of several dozens of additional LHC search card 



14 



files is provided along with the online program distri- 
bution package [18]. Users who create interesting new 
search cards are encouraged to submit them to the au- 
thor for inclusion with future releases. 



A. CMS Hadronic Multijet Search 



****** cut_card.dat 2.0 ****** 

* ASA : CMS Hadronic Multijets 

* CMS PAS SUS-11-003 k SUS-09-001 
*** Object Reconstruction **** 
0BJ_PH0 = PTM:25, CUT: [0,0] 
0BJ_ELE = PTM:10, CUT: [0,0] 
0BJ_MU0 = PTM:10, CUT: [0,0] 
0BJ_JET = PTM:30 

0BJ_JET_001 = SRC: +000, PTM:50, CUT: 9 
0BJ_JET_002 = SRC: +001, PTM:100, CUT: 2 
0BJ_JET_003 = SRC: +002, PRM : [0.0,2.5,1] , CUT:1 
0BJ_JET_004 = SRC: +001, PRM:3.0, CUT:[0,0] 
0BJ_JET_005 = SRC: +001, FEM:0.9, CUT: [0,0] 
****** Event Selection ******* 
EVT_MET_001 = JET: 000 
EVT_MET_002 = JET: 001, CUT: 100 
EVT_MHT_001 = JET: 001, CUT: 375 

EVT_RET_001 = NUM:002, DEN : 001 , CUT : [0 . , 1 . 25] 
****************************** 

V J 

Card A: CMS Hadronic Multijet Search [H Ojjj 

The first card file scenario to be documented here rep- 
resents a search by the CMS collaboration for SUSY sig- 
nals in purely Hadronic multijet events |19l . 1 351 ] . The 
procedure specified in the two CMS references differs in 
various minor details, and the provided card file thus 
represents something of a hybrid treatment. Following a 
point of some persistent interest to the author's research 
collaboration fl3L[l4l |. a high jet multiplicity threshold re- 
quirement is enforced by the current card that does not 
appear in the original CMS sources. This example like- 
wise omits a cut on the ai^ statistic of Eq. IT9|) , in keep- 
ing with the previously described difficulties facing appli- 
cation of that measure to extreme multijet content. The 
resulting selection is appropriate, for example, for com- 
parison against the data plotted in Figure lb of Ref. (3|| , 
which shows the event distribution per jet count, up to 
a maximum of 12, prior to application of the a T filter. 

The object reconstruction begins in line (5) by discard- 
ing events that possess photons meeting a 25 GeV min- 
imal threshold on the transverse momentum P T . Like- 
wise, events with either electrons or muons carrying a 
P T at or above 10 GeV are subsequently eliminated. The 
remaining object specifications are dedicated to the parti- 
tioning of hadronic jets. Line (8) prescribes a soft zeroth 
order cut of P T > 30 GeV on the transverse momen- 
tum of jets described in the input . Ihco event, which will 
be inherited by all subsequent groupings. A first posi- 
tively indexed partition, sourced from the prior zeroth 
order designation, is created from the jet object subset 
possessing a harder transverse momentum selection of 
P T > 50 GeV, and events are rejected that do not pos- 
sess at least 9 such jets. A second subdivision further 
requires the presence of at least two leading hard jets 



with P T > 100 GeV. The third jet reconstruction index 
extends this sequential processing of the object chain to 
require a pseudo-rapidity magnitude \rj\ < 2.5 for the 
leading jet; activation of the third position "stop" flag 
within the PRM parameter forces the clustering of objects 
to terminate upon the first selection failure (in descend- 
ing P T order), guaranteeing that a non-zero final count 
does indeed include the single hardest jet. The final two 
groupings in lines (12,13) are sourced retroactively from 
the previously defined P T > 50 GeV set in line (9), and 
further isolate jets with a pseudo-rapidity \rj\ > 3.0 or an 
electromagnetic fraction £ > 0.9, rejecting the event as a 
whole if any such surviving objects are identified. 

The event selection opens in lines (15,16) with miss- 
ing transverse energy definitions for two source popu- 
lations consisting, in turn, of all soft P T > 30 GeV 
jets and the harder P T > 50 GeV subset. The latter 
of these groupings, which is guaranteed from before to 
possess at least 9 objects with satisfactory rj and £ val- 
ues, must yield a missing energy of at least 100 GeV 
for the event to continue. Subsequently, the same set of 
harder jets is checked for its scalar transverse momen- 
tum sum iJ T , where a minimal value of 375 GeV is re- 
quired. Finally, an upper bound on the ratio of hard-to- 
soft jet missing energies (^ , ^, ard /^^? ft ) < 1.25 is imposed. 
If one wishes to reinstate the CMS ocl^ > 0.55 selec- 
tion, this is easily accomplished by addition of the line 
"EVT_ALT_001 = JET:001, CUT:0.55". Since the appropri- 
ate MET and MHT factors may be reconstructed from the 
provided jets they need not be supplied explicitly. 



B. ATLAS Hadronic Multijet Search 



1 ****** cut_card.dat 2.0 ****** 

2 * ATLAS Hadronic Multijets C7J80) 

3 * ATLAS-C0NF-2012-037 

4 *** Object Reconstruction **** 

5 0B J_ALL = PRM: [0.0,4.5] 

6 0BJ_ELE = PTM:20, PRM: [0.0,2. 47] 

7 OBJJTOO = PTM:10, PRM : [0.0,2.4] 

8 0BJ_JET = PTM:20, PRM: [0 . , 2 . 8] 

9 0BJ_LEP_001 = SRC:+000, EMT:+1 

10 0BJ_JET_002 = SRC:+000, CMP: +001, CDR:0.2 

11 0BJ_LEP_003 = SRC:+000, CMP:+002, EMT:-3, CDR:0.4, CUT: [0,0] 

12 0BJ_JET_003 = SRC: +002, PTM:40 

13 0BJ_JET_004 = SRC: +003, PTM:80, CUT: 7 

14 ****** Event Selection ******* 

15 EVT_MHT_001 = JET: 003 

16 EVT_RHT_001 = NUM:000, DEN: 001, CUT: 4.0 

17 ****************************** 



Card B: ATLAS Hadronic Multijet Search [13 

The second card file example to be presented repre- 
sents a parallel search by the ATLAS collaboration for 
SUSY signals in purely Hadronic multijet events [38j | . 
The specific selection strategy to be documented is re- 
ferred to as "7j80", in reference to the required jet count 
and transverse momentum per jet. Though targeting an 
essentially identical event profile to that just described, 
the mechanism of the adopted selections is distinct. 
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The object specification begins in line (5) with a global 
restriction on the pseudo-rapidity magnitude \rj\ < 4.5 
that will be enforced on all subsequently inherited ob- 
jects. The light lepton (electron, muon) populations are 
then filtered according to transverse momentum P T > 
(20, 10) GeV and pseudo-rapidity \t)\ < (2.47, 2.40). Like- 
wise, the zeroth jet classification is restricted to the ob- 
ject subset with P T > 20 GeV and \t]\ < 2.8. In line 
(9) a lepton reconstruction is sourced from the individual 
merged lepton flavors, and filtered specifically to con- 
tain only electrons. The inherited jet object set is then 
compared against the prior electron grouping for angu- 
lar proximity, and only those jets that are isolated by 
AR > 0.2 radians from the nearest electron are retained. 
Next, a fresh assembly of leptons is gathered in line (11) 
that specifically exempts tau flavor content. It is evalu- 
ated for angular proximity against all jets surviving the 
filter defined previously in line (10), and the event is dis- 
carded as a whole if any well isolated leptons, having 
AR > 0.4 with respect to all jets, persist. An interme- 
diate transverse momentum cut of P T > 40 GeV is then 
applied to the continued jets. Immediately following is a 
further hard classification of P T > 80 GeV, which must 
be survived by at least 7 jets for the event to proceed. 

The event selection phase initiated in line (15) is rather 
compact, firstly defining a scalar transverse momentum 
sum H^! ts for the intermediate jet classification, and 
secondly enforcing a relative significance boundary of 
^tVVI-^t* 8 ] — ^ GeV 1 / 2 on the missing transverse en- 
ergy $ T computed from all objects that were consistent 
with the primary pseudo-rapidity filter. 

C. CMS Dilepton and B-Tagged Jet Search 



****** cut_card.dat 2.0 ****** 

* CMS SSAF Dileptons and Bjets (SRI) 

* CMS PAS SUS-11-020 

*** Object Reconstruction **** 

0BJ_ELE = PRM: [1.566,1.422] 

0BJ_LEP = EMT:-3, PTM:10, PRM: [0.0,2.4] 

0BJ_JET = PTM:40, PRM: [0 . , 2 . 5] , CUT:2 

0BJ_LEP_001 = SRC:+000, SGN:+1, PTM : [UNDEF ,20] 

0BJ_LEP_002 = SRC: +000, SGN:-1, PTM : [UNDEF ,20] 

0BJ_LEP_003 = SRC: [+000,-001] 

0BJ_LEP_004 = SRC: [+000,-002] 

0BJ_LEP_005 = SRC: [+000,-001,-002] 

0BJ_JET_001 = SRC: +000, HFT : 1 , CUT: 2 

0BJ_DIL_001 = LEP:003, DLS:-1, DLF : 1 , DMI : [76 , 106] , CUT: [0,0] 

0BJ_DIL_002 = LEP:004, DLS:-1, DLF : 1 , DMI : [76 , 106] , CUT: [0,0] 

0BJ_DIL_003 = LEP:005, DLS:+1, DLF : , DMI:8, CUT:1 

****** Event Selection ******* 

EVT_MET = CUT: 30 

EVT_MHT_001 = JET: 000, CUT: 80 

****************************** 

V J 

Card C: CMS Dilepton and B-Tagged Jet Search [H| 

The next selection card deconstruction confronts an 
intricate CMS search strategy formulated around the in- 
variant mass of a same-sign, any-flavor dilepton pair, 
and requiring the presence of jets with heavy flavor tag- 
ging [39(. In particular, the example documents a region 



of the selection space referred to as "SRI". 

The object reconstruction begins in line (5) by dis- 
missing electrons that are located within an angu- 
lar dead region at the CMS barrel-endcap calorimeter 
seam corresponding to the pseudo-rapidity magnitude 
1.422 < \rj\ < 1.566 range. Next, there is a specifica- 
tion of the zeroth lepton object set that rejects intrusion 
by the tau, and places outer limits on the light flavor 
transverse momentum P T > 10 GeV and pseudo-rapidity 
\rj\ < 2.4. Likewise, the zeroth jet classification, for which 
surviving events must provide at least two matching ob- 
jects, is bounded by P T > 40 GeV and \r/\ < 2.5. In lines 
(8,9) a pair of signed lepton groupings is sourced from the 
established zeroth order set, one entirely positive and one 
entirely negative, each of which caps the transverse mo- 
mentum at P T < 20 GeV. Three additional lepton parti- 
tions are then recursively defined by the selective veto of 
candidates from the zeroth sourcing pool that also appear 
within one or either of the signed intermediate P T pairs; 
the result is a group of softer negative P T > 10 GeV lep- 
tons mixed with harder positive P T > 20 GeV leptons, 
a similar group with precisely the reversed charge versus 
transverse momentum relationship, and a group of mixed 
sign leptons that all respect the harder P T > 20 GeV clas- 
sification. A single indexed jet partition is reduced from 
the zeroth classification in line (13) by specifying a loose 
selector on the heavy flavor tag, and events without at 
least 2 such tagged jets are are cut. The final three ob- 
ject reconstructions are of the dilepton variety. The first 
two reject events containing an opposite-sign same-flavor 
dilepton configuration that yields an invariant mass M ' , 
as in Eq. ©, within 15 GeV of M z ~ 91 GeV, if at least 
one member of the pair is of the harder P T > 20 GeV ob- 
ject variety. The final object specification line requires 
the presence of at least 1 same-sign any-flavor dilepton 
pair with an invariant mass M 1 ' 1 > 8 GeV, where each 
lepton is above the individual P T > 20 GeV threshold. 

The event selection phase prescribed in lines (18,19) is 
quite straightforward, consisting of a cut > 30 GeV 
on the missing energy composed from all primary event 
objects, and a cut H^ ts > 80 GeV on the scalar trans- 
verse momentum sum for the zeroth assembly of jets. 

D. ATLAS Jets and Isolated Lepton Search 

The second ATLAS collaboration card file scenario 
to be documented represents a multi-level selection em- 
ployed to search for SUSY in final states with jets, miss- 
ing transverse momentum, and a single isolated lep- 
ton [2l|. An important role is played in this procedure 
by the transverse mass statistic of Eq. (fT"5l) . The exam- 
ple provided corresponds specifically to the 3-jet selection 
sub-strategy. Most ATLAS SUSY searches share a com- 
mon core of generic low-level object reconstruction filters, 
and the present card file specification does indeed bear 
several initial structural similarities to that featured in 
the previous ATLAS case study. 
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****** cut_card.dat 2.0 ****** 

* ATLAS Jets and Lepton (3J1L) 

* ATLAS-CONF-2012-041 

*** Object Reconstruction **** 
OBJ.ALL = PRM: [0.0,4.9] 
0BJ_ELE = PTM:10, PRM: [0 . , 2 . 47] 
0BJ_MU0 = PTM:10, PRM: [0 . , 2 . 4] 
0BJ_LEP_001 = SRC:+000, EMT:+1, PTM:25 
0BJ_LEP_002 = SRC : +000 , EMT:+2, PTM:20 

0BJ_JET_002 = SRC:+000, CMP:+001, PTM:20, PRM : [0.0,4.5] , CDR:0.2 
0BJ_LEP_003 = SRC: [+001, +002] , CMP:+002, CDR:0.4, CUT: [1,1] 
□BJ_JET_003 = SRC: +002, PTM:25, PRM : [0 . , 2 . 5] , CUT: 3 
0BJ_LEP_004 = SRC: [+000,-003] , EMT:-3, CUT : [0,0] 
0BJ_JET_004 = SRC:+003, CUT : [3 ,UNDEF , -1] 
□BJ_JET_005 = SRC: +003, PTM:80, CUT: [0,3] 
0BJ_JET_006 



SRC: +005 
****** Event Selection ******* 
EVT_MET = CUT: 250 
EVT_MHT_001 = LEP:003, JET: 004 
EVT_MEF_001 = MET: 000, MHT:001 
EVT_REF_001 = NUM:000, DEN : 001 , 
EVT_LTM_001 = LEP:003, MET: 000, 
EVT_MHT_002 = LEP:003, JET: 003 
EVT_MEF_002 = MET: 000, MHT:002, 



PTM:100, CUT:1 



CUT : . 3 
CUT: 100 



Card D: ATLAS Jets and Isolated Lepton Search |2l[ 

The object reconstruction begins in line (5) by enforc- 
ing a primary upper limit on the pseudo-rapidity mag- 
nitude \rj\ < 4.9 of all input objects. The light lepton 
(electron, muon) populations are then filtered according 
to transverse momentum P T > 10 GeV and pseudo- 
rapidity \rj\ < (2.47,2.40). In the next pair of lines a 
slightly harder variant of each of these two distinct lep- 
ton flavor populations is indexed relative to the implicit 
zeroth lepton composite, raising the lower limits on trans- 
verse momentum to 25 GeV and 20 GeV, respectively. 
An initial jet classification in line (10) is likewise sourced 
from its corresponding inclusive zeroth grouping, enforc- 
ing P T > 20 GeV and |ry| < 4.5 limits, and rejecting 
objects that are poorly isolated (AR < 0.2 radians) from 
members of the preceding electronic lepton partition de- 
fined in line (8). The electron and muon forks are sub- 
sequently rejoined into a unified object source for the 
next instruction line, which filters against a minimal iso- 
lation requirement of AR > 0.4 radians from elements of 
the just established jet reconstruction, and rejects events 
that do not afford 1, and only 1, compliant object. Af- 
ter this, a slightly harder and substantially more central 
jet partition is sourced from its direct numerical prede- 
cessor, with P T > 25 GeV and \r]\ < 2.5 limits, requir- 
ing a minimum of three such jets for event continuation. 
Line (13) marks a final return to the leptonic analysis, 
prescribing a cut on events where any non-tau flavored 
object from the zeroth classification proved too soft, for- 
ward or poorly isolated to pass through the remainder 
of the described selection cascade. The last three ob- 
ject reconstruction commands are additional jet assem- 
bly protocols. Picking up where the prior jet selection 
left off, the first of these lines imposes a redundant lower 
bound of 3 on the surviving jet content, but sets a value 
of for the third position "clip" input to the CUT pa- 
rameter, which prevents any jets beyond the minimally 
required count from appearing within the indexed object 
partition. The next instruction, which rejects events pos- 



sessing more than 3 hard jets with P T > 80 GeV, returns 
for its sourcing to the classification defined in line (12) 
rather than perpetuating the sequential filtering pattern; 
this discontinuity is essential, given that the immediately 
prior classification was capped at a 3 jet maximum. Re- 
suming a stepwise object flow, the final jet specification 
requires the event to contain at least 1 jet meeting the 
P T > 100 GeV threshold. 

The event selection phase opens in line (18) by plac- 
ing a lower bound on the inclusive missing energy 
250 GeV. Next, an indexed scalar transverse energy sum 
specification H^ ct is created from back-references to 
the single surviving lepton and the capped partition of 
3 jets, as defined in lines (11,14), respectively. Like- 
wise, an indexed effective mass specification is es- 
tablished as a sum of the prior two quantities. Subse- 
quently, a limit /{JPt + H^, 3)et ) > 0.3 is placed on 
the corresponding missing energy to effective mass ratio. 
Proceeding on, a cut > 100 GeV is instituted in 

line (22) on the transverse mass composed from the iso- 
lated lepton and the event missing momentum vector. 
In the next-to-final instruction, a more inclusive trans- 
verse energy sum H^) ets is defined that incorporates the 
lepton and all jets (not only the leading 3) within the 
P T > 25 GeV and \rj\ < 2.5 reconstruction from line (12). 
Finally, a lower bound is imposed on the effective mass 
(ff + P4 jots ) > 1,200 GeV formed by combining the 
previous scalar sum with the event missing energy. 



E. CMS Razor Variable Supersymmetry Search 



****** cut_card.dat 2.0 ****** 

* CMS Razor ELE Box CSR6) 

* CMS PAS SUS-12-005 
*** Object Reconstruction **** 
0BJ_ELE = PRM: [1.566,1.422] 
0BJ_MU0 = PRM: [0.0,2.4] 

0BJ_LEP = EMT:-3, PTM:10, PRM: [0.0,2.5] 
0BJ_JET = PTM:60, PRM : [0.0,3.0] 
0BJ_LEP_001 = SRC:+000, EMT:+1 
0BJ_LEP_002 = SRC: +000, EMT:+2 

SRC: +002, ETR: [0.00 
= SRC: +003, PTM:12, 
= SRC: +001, PTM:20, 
SRC: +003, CUT: [0,0] 
SRC:+002, PTM:15, CUT: [0,0] 
SRC: +002, CUT: [0,1], ANY: [006,007] 
= SRC: +001, PTM:20, CUT: [0,0] 
SRC:+001, CUT: [0,1], ANY:009 
SRC: +003, PTM:12 
SRC: +001, PTM:20 
****** Event Selection ******* 
EVT_JRM_001 = LEP:000, JET:000 
EVT_ALR_001 = LEP:000, JET: 000 



0BJ_LEP_003 = 

# 0BJ_LEP_004 

# 0BJ_LEP_005 
0BJ_LEP_006 = 
0BJ_LEP_007 = 
0BJ_LEP_008 = 

# 0BJ_LEP_009 
0BJ_LEP_010 = 
0BJ_LEP_011 = 
0BJ_LEP_012 = 



0.27] , PRM: [0.0,2.1] 
CUT : [0 , 0] 

CUT : [0 , 0] , ANY : 004 



CUT 
CUT 



[0,0] 



CUT: [450,1000] 
MET:000, CUT: [0.30,0 



50] 



Card E: CMS Razor Variable Supersymmetry Search I3ll, |32j 

The last card file scenario to be documented corre- 
sponds to a CMS search for SUSY using the razor vari- 
ables [U [HJ . Additional information on the active lep- 
ton reconstruction stra tegy and isolation requirements 
may be found in Ref. [20]. The selection strategy to 
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be highlighted hierarchically divides the harvested event 
content into one of several "boxes" according to the lep- 
ton (or dilepton) flavor content, or the lack thereof. To 
keep the various groupings disjoint, events satisfying the 
requirements of multiple boxes are uniquely sequestered 
upon their first match against an ordered sequence of 
clustered selection criteria. Conversely, the assignment 
of an event into one of the latter classification stages re- 
quires demonstrating that none of the prior available box 
qualifications were successfully matched. In general, each 
box will possess multiple positive attributes that must all 
be satisfied for inclusion of the event; therefore, a nega- 
tion of this inclusion is facilitated by falsifying any single 
member of the composite list of properties. The required 
logic may be implemented within the card file description 
by application of the ANY parameter to link related filters. 

The example provided corresponds specifically to sig- 
nal region 6 of the single isolated electron box from 
Ref. [32j. The object reconstruction begins in line (5) 
with the exclusion of electrons located within the pseudo- 
rapidity range corresponding to the barrel-endcap de- 
tector gap. Next, muons are required to be contained 
within \rj\ < 2.4, and a combined electron-muon lep- 
ton classification is composed with transverse momentum 
P T > 10 GeV and the slightly weaker pseudo-rapidity 
limit of \rj\ < 2.5. Similarly, the zeroth jet classifica- 
tion is defined with P x > 60 GeV and \r/\ < 3.0. In 
lines (9-11) the electron and muon objects are resepa- 
rated, and a tightened muon object grouping is instanti- 
ated with an explicit bound on the transverse momentum 
isolation ratio £ < 0.27 and a more central track orienta- 
tion \rj\ < 2.1. Lines (12,13) are the first content negation 
grouping to be tethered by the ANY keyword. The two 
members of this construct are actually made redundant 
by subsequent instructions, and, although their content is 
retained for pedagogical purposes, the leading comment 
markers # will suppress interpretation of the associated 
commands. In concert, events would be retained by this 
pair of instructions that lacked either a tight muon with 
P T > 12 GeV or a hard electron with P T > 20 GeV, 
both of which are mandated by the first box definition. 
The next three lines similarly retain events disqualified 
from inclusion in the second box, which is defined by the 
presence of at least one tight muon, at least one harder 
muon with P T > 15 GeV, and at least two muons over- 
all. Lines (17,18), one of which is again duplicative, would 
act collectively to exclude events that lack either a sin- 
gle electron with P T > 20 GeV or two electrons over- 
all. Line (19) constitutes a negation of the fourth lepton 
box criterion, namely the presence of at least one tight 
muon with P T > 12 GeV. Since this is a precise dupli- 
cation of line (12), which is logically coupled to line (13), 
there is no need to explicitly execute either of the prior 



tests. The final object reconstruction command grants 
otherwise successful events admission into the fifth lep- 
tonic box according to the presence of at least one hard 
electron with P T > 20 GeV. Surviving objects will thus 
necessarily fail the test embodied in line (17), which may 
consequently be omitted from its cluster. The same holds 
true for line (13), although its action was already blocked 
for independent reasons. 

The event selection phase prescribed in lines (22,23) of 
this card is quite simple, isolating a rectangular area of 
the razor plane with 450 GeV < < 1000 GeV and 
0.30 < ap J < 0.50 that is sourced from the inclusive miss- 
ing transverse energy and the full retained classifications 
of leptons and jets. 



V. AFTERWORD 

A new computer program named CutLHCO has been 
introduced, whose function is the implementation of 
generic data selection cuts on collider event specification 
files in the standardized .Ihco format. This software is 
intended to fill an open market niche for a lightweight 
yet flexible "consumer-level" alternative to the Root |10| 
data analysis framework. The primary envisioned appli- 
cation is as a filter on output produced by the PGS4 Q 
and Delphes @ detector simulations, which are them- 
selves lightweight alternatives to the Geant4 [7j based 
solutions favored by the large LHC experiments. All 
process control instructions are provided via a compact 
and powerful card file input syntax that efficiently facil- 
itates the reasonable approximation of most event selec- 
tion strategies and specialized discovery statistics com- 
monly employed by the CMS and ATLAS collaborations. 
The structure, function, invocation and usage of the most 
recent CutLHCO 2.0 program version has been docu- 
mented thoroughly, including a detailed deconstruction 
of several example card file specifications. The associated 
software is now available for free public download [18]. 
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Appendix: Card File Specification Synopsis 
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********* CutLHCO CARD SYNOPSIS Version 2.0 ************************************************************************* 

NAME = KEY. 1: VALUE. 1, KEY_2:VALUE_2; Multiple values go in square brackets; KEY : [VALUE_A , VALUE_B , VALUE_C] * 

Values default to zero-indexed option or UNDEF; Selectors with defined CUTs appear as output statistics * 

0BJ_ALL # Primary (first read access) classification for all objects; Unclassified objects discarded; (PTM, PRM) * 

0BJ_PH0 # Exclusive photon object classification; (PTM , PRM , CUT , JET) * 

0BJ_ELE # Exclusive solo electron object classification; Unclassified objects discarded; (SGN, PTM, PRM, CUT, JET) * 

0BJ_MU0 # Exclusive solo muon object classification; Unclassified objects discarded; (SGN, PTC, ETR, PTM, PRM, CUT, JET) * 

0BJ_TAU # Exclusive solo tau object classification; Unclassified objects discarded; (SGN, PTM, PRM, CUT, JET) * 

0BJ_LEP # Zeroth lepton object classification; Default input to latter stages; (EMT, SGN, PTC, ETR, PTM, PRM, CUT) * 

0BJ_JET # Zeroth jet object classification; Default input to latter stages; (HFT , FEM , TRK , MUO , PTM , PRM , CUT) * 

0BJ_LEP_N # Ordered lepton object classification; Index N > 0; (SRC, CMP, EMT, SGN, PTC, ETR, PTM, PRM, CDR, SDR, CUT, ANY) * 

0BJ_JET_N # Ordered jet object classification; Index N > 0; (SRC, CMP, HFT, FEM, TRK, MUO, PTM, PRM, CDR, SDR, CUT, ANY) * 

0BJ_DIL_N # Indexed dilepton object classification; Index N > 0; (LEP , DLS , DLF , DMI , CUT) * 

SRC : Source jet/lepton object classifications; [+/- Indices < N] ; >=0 : Includes ; <0:Excludes (Dominant) * 

CMP : Complementary lepton/jet object classifications; [+/- Indices < N] ; >=0 : Includes ; <0:Excludes (Dominant) * 

LEP : Local inclusive lepton object classification index; Index N >= * 

EMT : Accepted lepton flavors; 0:A11; +/-l:Elec; +/-2:Muon; +/-3:Tau; <0: Excludes * 

SGN : Accepted lepton sign; 0:A11; +l:Positive; -l:Negative * 

PTC : Accepted muon R=4 cone transverse momentum of adjacent objects range in GeV; [Min, Max] * 

ETR : Accepted muon adjacent transverse calorimeter energy & track momentum unitless ratio range; [Min,Max] * 

HFT : Accepted jet heavy flavor tagging; 0:A11; l:Loose (or Tight); 2:Tight (Only) * 

FEM : Accepted jet dimensionless electromagnetic fraction range; [Min.Max] * 

TRK : Accepted jet object composite track count range; [Min, Max] * 

MUO : Accepted jet integration of adjacent poorly isolated muons count range; [Min, Max] * 

PTM : Accepted transverse momentum magnitude range in GeV; [Min, Max, Sort] ; Sort=l : Ascending; Default : Descending * 

PRM : Accepted pseudo-rapidity modulus range in radians; [Min, Max, Stop] ; Stop=l : Terminate upon failure 

CDR : Accepted complementary inter-object Delta-R proximity range in radians; [Min, Max] 

SDR : Accepted source intra-object Delta-R proximity range in radians; [Min, Max] 

DLS : Accepted dilepton pair sign classification; -l:0pposite; 0:Any; +l:Same 

DLF : Accepted dilepton pair flavor classification; 0:Any; l:Same 

DMI : Accepted dilepton invariant mass range in GeV; [Min, Max, Sort] ; Sort=l : Ascending; Default : Descending 

CUT : Event cut if object count is outside given range; [Min, Max, Clip] ; Clip=-l:Clip to min; Clip=+l:Pad to max * 

ANY : Events passing any defined current or prior listed jet/lepton classification cut pass all; [Indices < N] * 

JET : Objects failing exclusive classification are recast as jets; l:True; Default :False * 

********* Event Selection ******************************************************************************************* 

EVT_CAL # Calorimeter based missing transverse energy event specification in GeV; Referenced as MET_N = -1; (CUT) * 

EVT_MET # Zeroth missing transverse energy (MET) event specification in GeV; Includes all objects; (CUT) * 

EVT_MHT # Zeroth scalar transverse momentum sum (MHT) event specification in GeV; Includes all objects; (MAS, CUT) * 

EVT_MEF # Zeroth effective mass (MEF = MET + MHT) event specification in GeV; Includes all objects; (CUT) * 

EVT_MET_N # Indexed MET event specification in GeV; Uses specified leptons & jets; Index N > 0; (LEP, JET, CUT) * 

EVT_MHT_N # Indexed MHT event specification in GeV; Uses specified leptons & jets; Index N > 0; (LEP, JET, MAS, CUT) * 

EVT_MEF_N # Indexed MEF event specification in GeV; Sum of specified MET and MHT; Index N > 0; (MET , MHT , CUT) * 

EVT_RET_N # Indexed unitless ratio (MET_NUM / MET_DEN) of specified factors; Index N > 0; (NUM,DEN,CUT) * 

EVT_RHT_N # Indexed ratio (MET_NUM / sqrt( MHT_DEN )) in GeV~{l/2} of specified factors; Index N > 0; (NUM.DEN.CUT) * 

EVT_REF_N # Indexed unitless ratio (MET_NUM / MEF_DEN) of specified factors; Index N > 0; (NUM,DEN,CUT) * 

EVT_DET_N # Indexed absolute vector difference (MET_TW0 - MET_0NE) of specified factors; Index N > 0; (ONE, TWO, CUT) * 

EVT_LTM_N # Indexed leading lepton & specified missing energy transverse mass in GeV; Index N > 0; (LEP, MET, CUT) * 

EVT_JSM_N # Indexed specified jet s-transverse mass MT2 in GeV; Index N > 0; (JET, CUT) * 

EVT_JZB_N # Indexed leading dilepton fc specified jet Z-balance in GeV; Index N > 0; (JET,DIL,CUT) * 

EVT_JRM_N # Indexed specified lepton & jet razor mass M_R in GeV; Index N > 0; (LEP, JET, CUT) * 

EVT_ALR_N # Indexed specified lepton, jet & optional MET unitless alpha_R ratio; Index N > 0; (LEP, JET, MET, CUT) * 

EVT_ALT_N # Indexed specified jet & optional MET, MHT unitless alpha_T ratio; Index N > 0; (JET, MET, MHT, MAS, CUT) * 

EVT_MDP_N # Indexed specified jet & optional MET minimal delta phi angle in radians; Index N > 0; (JET, MET, CUT) * 

EVT_BDP_N # Indexed specified jet & optional MET biased delta phi angle in radians; Index N > 0; (JET, MET, CUT) * 

LEP : Local inclusive lepton object classification index; Index N >= * 

JET : Local inclusive jet object classification index; Index N >= * 

DIL : Local inclusive dilepton object classification index; Index N >= * 

MET : Missing transverse energy specification index; Index N >= -1; Defaults to local objects when optional * 

MHT : Scalar transverse momentum sum specification index; Index N >= 0; Defaults to local objects when optional * 

MEF : Effective mass event specification index; Index N >= 0; Defaults to local objects when optional * 

MAS : Masses of individual objects are retained during intermediate computations; l:True; Default : False * 

NUM : Numerator content specification index; Refers to one of (MET, MHT, MEF) by context; Index N >= -1 or * 

DEN : Denominator content specification index; Refers to one of (MET, MHT, MEF) by context; Index N >= -1 or * 

ONE : First input content specification index; Refers to one of (MET, MHT, MEF) by context; Index N >= -1 or * 

TWO : Second input content specification index; Refers to one of (MET, MHT, MEF) by context; Index N >= -1 or * 

CUT : Event cut if associated value is outside given range in contextually established units; [Min, Max] * 



